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Abstract 

We construct general models for holographic superconductivity parametrized by three couplings 
which are functions of a real scalar field and show that under general assumptions they describe super- 
conducting phase transitions. While some features are universal and model independent, important 
aspects of the quantum critical behavior strongly depend on the choice of couplings, such eis the order 
of the phase transition and critical exponents of second-order phase transitions. In particular, we 
study a one-parameter model where the phase transition changes from second to first order above 
some critical value of the parameter and a model with tunable critical exponents. 
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1 Introduction 

Recently, the AdS/CFT correspondence [3 121 12 was applied to study strongly coupled systems which 
undergo a superconducting phase transition below a critical temperature [H [51 [SJ [71 |H] ■ On the field the- 
ory side, in these models superconductivity is characterized by the condensation of a composite charged 
operator for low temperatures T < Tc- In the dual gravitational description, the superconducting phase 
transition is represented by a transition from the black hole in anti-de Sitter space, with Hawking tem- 
perature equal to T, to a new solution with scalar "hair" , which is thermodynamically preferred below 
the critical temperature Tc- 

The models studied here are a generalization of the models introduced in [51 [TU]. In particular cases, 
they reduce to the 3-1-1 dimensional model studied in [TT], or the 4-1-1 dimensional model of |T2]. Certain 
properties of superconductors, like the London equation or infinite DC conductivity, follow from U{1) 
spontaneous symmetry breaking and are therefore universal (see e.g. |13j). Other features, like the 
number of vacua and the detailed dynamics of the phase transition (including critical exponents), may 
be model dependent, as will be illustrated by means of some examples. Therefore the family of models 
constructed here could be used as a set-up for model building. 

Another motivation in our construction is to have a framework where one can continuously interpolate 
between different models appeared in the literature. Finally, we also hope that a special choice of 
couplings could lead to simpler (perhaps analytical) solutions, and one example (albeit irrelevant for 
superconductivity) of a model where the equations of motions are reduced to first-order differential 
equations is given in appendix A. 

This paper is organized as follows. In section 2.1 we introduce the family of models and describe 
some basic conditions that the couplings should satisfy in order to have AdS/CFT duality at work and 
a consistent model of holographic superconductivity. In section 2.2 we write down the ansatz and the 
equations of motion. We also describe the boundary conditions that give rise to the hairy black hole 
solutions that lead to a U{1) spontaneous symmetry breaking in the boundary field theory. In section 
3 we present a preliminary discussion on how some aspects of the phase transitions should be encoded 
in the couplings. This is summarized by four conjectures. We then present the numerical results for the 
integration of the differential equations (including the back reaction) in three different models. These 
novel models are designed to test the conjectures made at the beginning of the section. The models of 
sect 3.1 and 3.2 give rise to features which already appeared in inequivalent models in [SllTO], like the 
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passage between first and second order phase transitions, and tunable critical exponents. The model 
of section 3.3 describes a physics similar to the model of Hartnoll, Herzog and Horowitz [7] (hereafter 
HHH model). Section 4 contains concluding remarks. Appendix A gives the equations in a slightly more 
general gauge and the expression for the conserved charge which is present for every model in our family. 
In this appendix we also show that for a particular choice of couplings there is an extra conserved charge 
which represents an integral of a combination of the equations of motion. This permits to write the 
equations of motion as first-order differential equations. Finally, appendix B contains a simple derivation 
of the AdS Reissner-Nordstrom (RNAdS) solution by the method of the superpotential. 



2 Generalized superconducting model 
2.1 The model 

Consider the following U{1) invariant d + 1 dimensional Lagrangian 



-gC^^l F{r^) R - iG(7?) F^^'-F^, - V{ri) - \H{^){diif - \jmd^e - A^f ] . (2.1) 



with /i, V ~ Q^l, ...,d. It represents the most general covariant Lagrangian with local U{\) symmetry, 
containing a metric g^,y, a U{\) gauge field A^, a real scalar field 77 and a Stiickelberg field 9, and terms 
with no more than two derivatives. The gauge transformations are the standard ones given by 

A,,^ A,,+d^,K , e^e + K, (2.2) 

with g^n and 77 being invariant. 

For d > 1, by a suitable Weyl rescahng 5^1, — > f{r])g^„ and a field redefinition r] — > K{r]), we can 
remove two of the couplings and put the Lagrangian in the following form 



^C^V^( R-\g{v) F^'^F,, + ^^^C/(r7) - \{diif \j{^)A,A^ ) , (2.3) 



where we have also fixed the gauge = and relabeled the couplings. The Lagrangian (|2.3p is our 
starting point for the construction of interesting models for holographic superconductors. Our aim is to 
understand the effect that G(r]), U{r]) and J^r]) have on physical quantities and on the dynamics of the 
phase transition. 

The requirement that the model exhibits holographic superconductivity implies a number of con- 
straints on the behavior of the couplings G, U, J. In particular, the couplings must be such that the 
theory contains the usual AdS solution. This requires that there is an extrema of U{ri) at which C/(r/) 
has a finite, positive value. By a suitable shift in r/, with no loss of generality this extrema can be placed 
at 7y = and U (0) can be set to 1 (any other value of U{0) > can be absorbed into a redefinition of the 
parameter L). Hence, at small 77, one must have the behavior U{ri) = 1 + 0{Tf). The potential could also 
be such that it admits other extrema. An example is the model with quartic potential U (rj) investigated 
in [5J[T3]. On the other hand, G{r]) and J{ri) must be positive to ensure that the kinetic terms for 
and 9 have the correct signs. 

We also request the theory to contain charged black hole solutions with AdS asymptotics and 77 = 
(here we will not discuss models with dilatonic black holes that may also exhibit interesting critical 
phenomena, see [19j for some examples). This implies that G(0) must be finite and non-zero and with 
no loss of generality we can take G(0) = 1 (which fixes the canonical normalization for A^). In addition, 
the theory must also contain black hole solutions with non-trivial scalar hair 77 with appropriate behavior 
at infinity, which implies that it must goes to zero in such a way that the dual field operator takes an 
expectation value (this point will be expanded below). In addition, J(ry) must vanish as r/ — )■ since we 
want U{1) to be unbroken at 7; = 0. The condition that the theory contains charged black hole solutions 
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with 77 = also implies that 9,j J and dr/G must vanish as 77 —> 0, in order to satisfy the 77 equation of 
motion. Thus, in this paper, the couplings will be required to have the following behavior at small 77: 

G(77) ^ l + crf + 0{rf), 

2 

Uiv) - ^-^^^V^+Oi,^), (2.4) 

J(77) - ^ o(^3) 

An important subclass are models with Z2 symmetry 77 —7/, for which the required behavior is 

G(77) = 1 + C772 + 0(77^) , 

? 

Uiv) - ^-^;^)V^+Oi,^), (2.5) 
Jiv) = qW + 0{v^). 

i.e. the expansion contains only even powers of ?/. The mass parameter rri^ is required to satisfy the 
Breitenlohner-Freedman bound for stability, tti^ > /A. It should be noted that it is not necessary to 
have g 7^ in order to have holographic superconductivity. In particular, in |2] it was shown that the 
choice J ~ ri'^ leads to superconductivity as well. 

Different choices of G , U, , J reproduce the various models that appeared in the literature. In 
particular, the HHII model [7 corresponds to taking d — 3 and the choice 

G(77) = 1 , U{v) = l + ^Tj\ J{rj) = q^rf . (2.6) 

For a general model with the asymptotic (|2.4p and tti^ — —2, the onset of instabilities will occur at the 
same critical temperature as in [7], since in this case the linearized Lagrangian is the same. 

Other c? = 3 models that appeared in the literature are obtained by the following identifications. The 
model of Franco et al [9] is the choice 

G(77) = 1 , C/(77) = 1 + 1 772 , Jin) = arbitrary . (2.7) 

The M-theory model of Gauntlett et al corresponds to taking 

G(77) = 1 , U{v)^^ cosh^ I (7 - cosh 77) , J(77) = ^ sinh^ 77 . (2.8) 

The string theory model of Gubser et al ^2j is obtained by considering the same Lagrangian (|2.3p with 
c? = 4 and setting 

G(77) = 1 , 1/(77) = i cosh^ ^ (5 - cosh 77) , JM = ^ sinh^ 77 . (2.9) 
2 2 

An even wider class of models can be obtained if one defines 77 = e"^, with the fundamental degree 
of freedom being ip, a real scalar field. In other words, if we assume that 7; > 0. This permits a more 
general dependence on 77, that includes non-integer powers. In particular, the couplings would now be 
required to have the asymptotic behavior 



G{r)) = 1 + 0(7/ 
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Uiv) - l-^^7(^^^ + 0(77-)> (2.10) 
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with ai, a2, being any real numbers subject to the condition ai > l,a2,3 > 2 (the condition ai > 1 
comes from the requirement that the RNAdS black hole with 77 = solves the 77 equation). A particular 
example of this is the d = 3 model investigated in [10] , which in our notation reads 

G{tj) = 1 , U{v) = 1 + i , J{r^) = rf + Ca^r^"^ + C4rj^ , 3 < a < 4 . (2.11) 

D 

This model exhibits interesting features like tunable critical exponents, i.e. which change by varying the 
parameter a. 

Thus the Lagrangian (|2.3p provides a universal structure that embodies various models studied in the 
literature and permits a more general approach to the phenomenology of superconductors. 



2.2 Holographic superconductivity 

We now set d = 3 and consider the following ansatz: 

ds^ = -g{r)e-^^''Ut^ + — +r^(dx^ + dy^), A ^ (t>(r)dt , V = v{r) ■ (2.12) 



The effective Lagrangian takes the form 

-^C = -2e-i{rg)' + -G{v)eic^'' + ^e-${/(r;) - -e'^gv'' + -eij{rj)ct>^ (2.13) 
and the equations of motion reduce to 

X' + y + T^e>^JivW = , (2.14) 



1 ^,2 , G{f]) ^^^,2 , .g' , 1 3 ^^^^^ , 1 „y;,^„^^2 

4 " 



7 V + ^ +ZZ + -~Tk Uirj) + — e^Jir,)^'' = , (2.15) 



G(ry) 






1 3 




rg 


J.2 ^2g 




x' 






(^ 


2 


G 


y 5G(77) 




x' 

2 


9 J 


2g ^ 



= 0, (2.16) 

^J?!)2=0. (2.17) 



In the appendix A we give the equations in a more general coordinate system and discuss the presence 
of Noether charges. 

The basic strategy to solve these equations is explained in detail in [71 18] . Here we summarize the 
main features. Black holes with regular event horizon are needed in order to put the field theory at finite 
temperature. We require that g has a simple zero at the event horizon r — r+, that is, g{r+) — but 
g'{r+) 7^ and finite. The Hawking temperature of the black hole is, 

THawk = ^9'{r)e-^'^-^/' . (2.18) 

On the other hand, g'{r^) is not an independent parameter in our gauge but it is determined by a 
combination of eqs. (|2.14l) and (12. 15^ . 

/ x\' 3r2 x^^, , r^G(ri) x ,,2 

[rge--j = —e--U{v) ^ e^0'' . (2.19) 

Substituting this value of g' into the Hawking temperature (I2.18P we find 

THawk = [e-^12U(v+) - L'e^G{v+) El) , (2.20) 
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where we have defined 



b'{r+) = E+, x{r+)=X+ 



(2.21) 



Solving the full set of equations involves numerical integration from the horizon to infinity. Therefore, 
in order to have a well-posed Cauchy problem, we have to specify the remaining initial data at r = r_|_. 
By expanding (|2.17p in the vicinity of the horizon one determines the value of ri'{r^) = rj'^ through the 
relation, 

ry;.g'(r+) + 1 e^^d,G{rj+) El + ^d^U{r,+ ) = . (2.22) 

In addition, the condition (^(r^.) = is needed to ensure that the gauge one- form is well-defined at the 
horizon (see [4] for a discussion). The solution is therefore specified by the values r+, 77+, Ej^ and x+- 
This set of parameters can be further reduced by using two scaling symmetries of the metric, the gauge 
field and the equations of motion. The scaling dimensions are summarized in the following table: 



Symm. 




t 


X 


r 


9 





I 


2 


1 











-1 


II 





-1 


-1 


1 


2 


1 



(2.23) 



The symmetry II can be used to set r_|_ = 1. The symmetry I can be used to set Xo = Oj where xo is tbe 
asymptotic value of x- Thus the initial value data is completely characterized in terms of two parameters, 
77+ and 

Because of the conditions (|2.4p . the asymptotic behavior of the solution will be essentially the same 
as in the HHII model. At infinity, the solution approaches the AdS geometry with 



r 

I2 



Using (|2.16p - p.l7p one can then find the asymptotic expressions for (f) and 77, 



r]{r) 



.(1) 



n(2) 



(2.24) 



(2.25) 



Here A± are the roots of A(A — 3) = w?. The gauge-gravity correspondence proceeds as usual (see e.g. 
[15l[T6]). The scalar field is dual to an operator Oa of dimension A equal to the larger root A+, while /i 
and p are interpreted as the chemical potential and the charge density of the dual field theory. In order 
to have spontaneous symmetry breaking in the field theory, we have to consider solutions with 77^^^ = 0. 
This condition is, in fact, equivalent to demanding that 0\, is not sourced. For m? in the interval 



and 0\ 



as both modes with t;^^) or 77 



(2) 



— 9/4 < 7)7,'' < —5/4 it is also consistent to swap the role of Oa 
are normalizable. 

The set of pairs (7;+, i?+) solving the condition r]'^^\r]+, E+) = define the condensed phase in which 
(Oa+) 0- As a result of this condition, the hairy black hole solution is fully characterized by a single 
parameter, which can be taken to be the temperature. 



The non-condensed phase is obviously defined by (Oa) — 0, i.e. 77 



(1) 



V 



(2) 



0. This condition 



implies 77 = 0, therefore this phase is dual to the usual AdS Reissner-Nordstrom solution 



1 



g = r 



P 

Ar4 



P 

4^2 



= P 



1 

r+ 



(2.26) 



in units where L — I (see appendix B for a derivation). 

Superconductivity only exists if the condensed phase is thermodynamically preferred, namely if the 
difference of free energy, 

A/(T) = /co„d(r) - /uncond(T) , (2.27) 

is negative for some range of temperatures. In the next section we will see typical examples where, at 
sufficiently low temperatures. A/ is negative and the system is in a superconductor state; then it increases 
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as the temperature is increased up to some critical value T = Tc at which A/(Tc) — 0. For temperatures 
above the uncondensed solution has less free energy and becomes the thermodynamically favorable; 
the system turns into a normal state. The behavior of A/ at T = Tc determines the type of phase 
transition, according to the standard classification. For our system the free energy can be expressed 
in terms of the energy density, chemical potential and density charge. The energy density of the field 
theory configuration is associated with the mass e of the black hole solution, which can be read from the 
asymptotic expression for the 00 component of the metric, 

e-X5 = r2-± + ... (2.28) 

A similar calculation as in [5] gives 

for all (G, J7, J) models. We will study the temperature dependence 
by working at fixed charge density p. 

3 Numerical Results for three examples 

General models with couplings that have the asymptotic behavior (j2.4p are naively expected to exhibit 
similar dynamics near the critical point, where t] is small and higher order terms in the expansion in 
powers of t] are in principle irrelevant. In particular, this expectation is partly confirmed by results 
obtained in [11015], which show that the details of the potential I] are irrelevant as far as the existence 
of a condensed phase is concerned. The mass parameter in [/ = 1 — m?rj^ /\2 + ... affects the critical 
temperature but the qualitative features of the transition are essentially the same |17| . Moreover, it was 
shown in [18 that for the HHH model the critical exponents take the standard mean field values. 

Nonetheless, this picture does not apply for all models. First, as shown in [S] for the J — rf model, 
for certain models the transition can be first order and occur at a finite (non-small) value of 7^0 The first 
example discussed here will illustrate this fact. It consists of a model obtained by a choice of non-trivial, 
one-parameter coupling G{ri) in the Maxwell kinetic term, while having the same couplings U{ri) and J{ri) 
as in the HHII model. The model exhibits a number of interesting features that we discuss. Second, there 
are also models exhibiting second-order transitions where critical exponents are modified. An example 
is the model (|2.1ip studied in [TU]. The critical exponent /3 for the order parameter (O2) ~ (7c — T)^ is 
given by the relation /3 = {a — 2)~^ . Although rj is very small near the transition, the term 77" in (|2.1ip 
is not sufficiently suppressed and it affects the critical behavior near Tc. 

By an intuition based on the few known examples, one may conjecture that 

1) For Z2-symmetric models with the asymptotic behavior (|2.5p . with q 0, there is a certain range 
of parameters where they undergo a second-order phase transition. They have the standard critical 
exponents predicted by mean field theory. 

2) In these models the critical temperature of a second-order phase transition is sensitive to the values 
of the parameters c, m, q and insensitive to the coefficients of the following powers. 

3) For models which are not Z2 symmetric, with the asymptotic behavior (j2.4l) . and <? 7^ 0, the critical 
exponents of second order transitions can be affected by the powers 0(?7'^) in (|2.4p but they are insensitive 
to the next terms starting with 0(77^). 

4) Models with ?7 > 0, i.e. having a real scalar field as a fundamental variable, defined by 77 = e''' > 0, 
can have the asymptotic behavior (I2.10p . In a regime of parameters where the couplings G, U and J 

^ In this paper, back reaction on the geometry is neglected -an assumption that can be justified to some extent by taking 
a suitably limit involving large C3 in a J = csr}^ model. Thus a natural question is whether these important changes in the 
dynamics of the phase transition might be an artifact of the no back-reaction approximation. We have checked that the 
study of the phase transition including the complete back-reaction reproduces essentially the same features observed in [9] , 
even for low values of C3 . 



(2.29) 

of the free energy in different models 



7 



are congruent to (|2.6p , these models have second-order phase transitions with critical exponents that can 
depend on the values of the parameters. 

The examples given below, with three different choices of G{ri), partly confirm this picture (and 
complement the study of [9l[TTl[T0], performed for G = 1 and different U{ri), J{ri) couplings). 

3.1 Model with G{r]) = (1 + br]'^y^ 
Consider a model with the following couplings 

The special value b = corresponds to the HHH model. The theory contains a variety of black hole 
solutions. The uncondensed 77 = phase is described as usual by the RNAdS solution (|2.26p . When the 
temperature is below some critical value Tc, one or more hairy black hole solutions appear, depending on 
the value of the parameter b. For b below some critical value 6cr ^ 3, there is a single black hole with 
hair. When b is greater than 6cr, new hairy black holes appear. Remarkably, the picture turns out to 
be qualitatively similar as in the model discussed in [10, with G = 1, U ~ 1 + ^ rj^, and J = + C477'*, 
with C4 playing the role of b. One difference is that in the present case the critical temperature depends 
on b (in consistency with conjecture 2), whereas for the model of fTU] with J — rj^ + C4,rf the critical 
temperature remains fixed and independent of C4 (which is also a prediction of conjecture 2). 



(a) (b) 

Figure 1: (a) Value of the condensate as a function of the temperature for the model with G = (1 + 677^)"^, 
U — 1 + J = Qif, for b = 0.5, 3, 5, 9 (from bottom to top, the b = 0.5 case being described by the 
undermost curve at T = 0). (b) Normalized free energy for the same model and same values of b (the 
b = 0.5 case being described by the rightmost curve). The lower branches of the 6 = 5 and & = 9 curves 
correspond to the upper branches in fig. la. 

Figure la shows the condensate (O2) as function of the temperature for qL ~ 3 and b — 0.5, 3, 5, 9. 
One can see that for b above a critical value &cr ~ 3 there is a new branch of the solution. Figure 2b 
shows the difference between the free energy of the solution and the free energy of the RNAdS black hole. 
From this figure one can observe a number of important features. 

For small values of &, the curves reproduce, as expected, the same qualitative features of the HHH 
model: there is only one hairy black hole solution, and the normalized free energy (i.e. the free energy 
of the hairy solution minus the free energy of the RNAdS solution) becomes negative exactly at T < Tc, 
showing that the black hole with scalar hair dominates the thermodynamics in this regime. The first 
derivative of the free energy is continuous at the transition, showing that the transition is second order. 
The condensate approaches zero as (O2) ^ {Tc — T)", with a~ 1/2 for all b < b^r- 

For b greater than &cr, the curves become multivalued, corresponding to the presence of two black 
hole solutions with different values of (O2) at the same temperature. Figure 2b shows that the upper 
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branch (i.e. with higher (O2)) has less free energy than the lower branch, and hence dominates the 
thermodynamics. 

The critical temperature Tc is the temperature at which the normalized free energy changes sign. 
Above Tc, the normalized free energy becomes positive and the RNAdS solution becomes thermodynam- 
ically favorable. As can be seen from figs, la, lb, this temperature is below the temperature at which 
the lower and upper branches of the solution join. This implies that for T > Tc the thermodynamically 
relevant solution is RNAdS, and the vacuum expectation value has a discontinuous jump from a non-zero 
positive value to (O2) = 0. 

For b > 6cr the normalized free energy is continuous but not differentiable at T = Tc; proving that 
the transition becomes first order. 



3.2 Model with €{7]) = (1 + 

This example illustrates a model with r/ > and represents another check of conjecture 4 (the first check 
is the model (|2.11l) of [10] with /?=(« — 2)^^). We have used the value b = 0.5, with three diff'erent 
values of a, a = 3.75, 3.5, 3.25. For simplicity -and to achieve greater numerical precision- the figures 
have been obtained ignoring back reaction, but we have verified that the incorporation of back reaction 
gives essentially the same results. For these values of the parameters the model has second-order phase 
transitions. Figure 2a shows the order parameter (O2) in terms of the temperature T for these three 
different values of a. We can see that the critical temperature is independent of a, Tc = 0.11843. Figure 
2b shows log(02) vs. log(l — T/Tc). Near the critical temperature we have the behavior 

(O2) - {Tc - Tf , (3.2) 

with [3 = 0.58, 0.67, 0.81 for a = 3.75, 3.5, 3.25 respectively. Remarkably, f3 is related to a by the same 
law as in the model (|2.1ip : 

P = {a-2)-^. (3.3) 

Despite this coincidence, the two models are inequivalent (though they seem to belong to the same 
universality class). As a check of conjecture 3, we have also studied the case a = 3, which gives /? = 1. 




T i-T/T, 
(a) (b) 



Figure 2: (a) Value of the condensate as a function of the temperature for the model with G = (1 + 
0.5??")"^, U = 1 + J = T]'^, for (from top to bottom) a = 3.75 (orange), a = 3.5 (red), a — 3.25 
(blue), (b) (O2) vs. (1 — T/Tc) in logarithmic scale, showing that different values of a have different 
slopes (conventions as in (a)). 
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3.3 Model with G{r]) = cosh 77 



Our third example illustrates the fact that there are as well non-trivial choices of G{ri) defining models 
which share essentially the same physics as the HHH model, even when the coupling becomes exponentially 
large at large 77. We consider the model defined by the couplings 

G = coshr], U = l + —, J^q^jf. (3.4) 

6 

In a sense, this model is complementary to the model of section 3.1, where G was suppressed at large 77. 
If r] is small in the whole space r+ < r < 00, then cosh 77 ^ 1 and the couplings are essentially the same 
as in the HHH model. If 77_|_ is large, then G becomes exponentially large near the horizon region and in 
principle one may expect some different physics. The results we have found, summarized by figures 3a 
and 3b, are as follows. 

The uncondensed 77 = phase is as usual described by the RNAdS solution. Below some critical 
temperature Tc ~ 0.188 we again find a hairy black hole solution. This critical temperature is slightly 
above the one of the HHH model. In general, we have checked that for a model with G{'q) = coshaTy the 
critical temperature depends on the coefficient a. This is consistent with the above conjecture 2. 

Figure 3a shows the value of the condensate as a function of the temperature, whereas figure 3b shows 
the normalized free energy. Below the critical temperature Tc, a black hole solution with scalar hair 
appears. Figure 3b shows that this has less free energy than the RNAdS solution, implying the existence 
of a phase transition. The free energy has a continuous derivative at the critical point, which shows that 
the transition is second order. A fit of the curve near the critical point shows that the critical exponent 
is given by the mean field value, (3 « 1/2, in consistency with conjecture 1. 




0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.08 0.10 0.12 0.14 0.16 0.18 0.20 

T T 

(a) (b) 

Figure 3: (a) Value of the condensate as a function of the temperature for the model with G — cosh?], 
[7=1 + ^, J = Or/^ . (b) Normalized free energy for the same model. 



4 Conductivity 

In this section we briefiy discuss the set-up for studying conductivity in the general model (|2.3p . We will 
consider d — 3 and G,U, J with the asymptotic conditions p.4|) with = — 2. Following ^ E] [2, we 
consider time-dependent perturbations of = ax{r) e~*'^* and gtx = fir) e^*'^*. In our general model, 
these fluctuations are governed by the following equations 

/'--/ + G <P'ax - . (4.2) 
r 
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Substituting the second into the first equation, we find 




(4.3) 



The asymptotic behavior of the perturbations is found to be 



a. = 40) + ^ + ... , f^r^f(o) + ^ + ... (4.4) 



r r 
The conductivity can then be found by the formula 



E (0) ' 



(4.5) 



where in the second equahty we have used the AdS/CFT dictionary. 

At this point, we can now examine under which conditions the general class of models (j2.3p lead to 
infinite DC (i.e. at cj = 0) conductivity. The emergence of infinite DC conductivity is a general feature 
that follows after the U{1) symmetry breaking due to a mass term for the photon. In the present case, 
the Maxwell equation has a term which is the one to be interpreted as a (London) current. As 

long as this is non- vanishing, one expects as usual infinite DC conductivity. The other terms in the above 
Maxwell equation are due to back reaction or interaction of the Maxwell field with the geometry. To be 
more precise, we need that a'x^ /o'x^ is a real, non- vanishing number in the limit that — 0, so that the 
imaginary part of the conductivity contains a pole, Im((T) ^ 1/lo. Then, by virtue of the Kramers-Kronig 
relation 

. / / XX 1 ^ /■ Re(cr(w'))dw' 

Im(a(a.)) ^ — V ^ ; , (4.6) 

one finds that a pole in Im(cr) is associated with a Dirac delta function in Re((T(w)). Typically, one 
expects that any J which vanishes at 77 = (i.e. in the normal, AdS Reissner-Nordstrom phase) and 
does not vanish for 77 ^ should lead to superconductivity, if it also leads to a phase transition like those 
described in the previous sections. 



5 Concluding remarks 

Summarizing, we have presented a large class of models which can be used as a phenomenological set-up 
for superconductivity or for model building. The general models interpolate between different models 
appeared in the literature. We have discussed general conditions that the couplings must satisfy for 
the theory to holographically describe superconducting phase transitions. In section 3 we have also 
conjectured how some aspects of the dynamics of the phase transition are encoded in the couplings 
G, J7, J. The concrete examples investigated here partly confirm this picture, but a more thorough 
investigation would clearly be necessary to substantiate these conjectures, possibly by analytic, rather 
than numerical, techniques. We have explored models with non-trivial G{rf). From a more fundamental 
standpoint, such models might arise from M-theory or string theory compactifications, or by quantum 
effects, since such couplings are compatible with U{1) local symmetry and general covariance. 

There are, of course, other aspects of superconducting phase transitions which have not been discussed 
here. There may be many surprises in the space of models parametrized by the couplings G(ry), U{ri), J{ri) 
and our analysis only explores a small corner of this space. In particular, one feature that emerges in 
certain regions of the parameter space for some models is the appearance of critical curves where the sign 
of the second derivative {O2) / dT"^ changes at some T < Tc- In some cases this (probably unphysical) 
feature is cured by the incorporation of back-reaction, but we have found cases where this feature persists 
even after including back reaction effects. 
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Another mteresting problem is to look for models with special couplings G{ri), U{rf), J{r]) which 
could allow one to find analytical solutions. The general system (|2.14P " (|2.17p of four coupled differential 
equations contains two first order equations and two second order equations (12.151) . Using the conserved 
charge discussed in the appendix A, one can replace the equation for (|2.16p by a first-order one. A 
natural question is then if there is some choice of couplings where all equations are first order. This 
would render the numerical analysis simpler and perhaps lead to an analytic treatment of at least some 
aspects of the phase transitions. A simple example is presented in appendix A, where we show that 
indeed there is one model which has an extra conserved charge and the equation for can be substituted 
by a first-order equation. Although this particular model does not seem to be relevant for holographic 
superconductivity, it would be interesting to see if there are other choices of 0(77), [/(yy), J(ri) where the 
equations can also be integrated to a first-order form (see also the superpotential method described in 
appendix B). 
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A General ansatz and conserved charges 

Consider the following ansatz 

ds^ = -g(r)e-^^'''>df + ^+ a^(r) (dx^ + dy^) , A = (j)(r)dt , 77 = r](r) . (A.l) 

which is slightly more general than the ansatz (|2.12p adopted in section 2.2. The Lagrangian (|2.3I) takes 
the form 

6a^ _x^,, , _x ,2 X 



V^C = 2e-ia'[{agy - agx') + yG(r7)e40'' + ^e~iu{v) ^e~igrj" + ^g^^Av)^' (A.2) 
The equations of motion are given by 

2a" + a'x' + + = , (A.3) 



1 ,2 G(t]) ^ ,2 a'q' a'^ 2a" 3 , , 1 ^ , . 9 . 
4 45 ag a a L^g ig'^ 



, /2a' x' , d,Gv'\ Jiv) 



2 G J gG{rj) 



= , (A.5) 



r," + v' [ — + e^^^G 0'' + ^d^U+^ e^d^J = . (A.6) 



a 2 gj • 2g ^"^ ^ L^g ^"^ 2g^ 
One can show, generalizing the similar derivation given in [14 , that the following quantity 



Qi = a2e>^/2 (a^ (^^) ' _ G^A (A.7) 



is a constant. This is a conserved charge associated with the following scaling symmetry 



a 



^c>, a^-, x^X-41og(c), (A.8 
c 
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5^5, V^V ■ (A.9) 

A natural question is if there is a particular choice of the couplings G, U, J for which the system has some 
extra conserved charge. We found that the system with 

6 = 6"", J = q^e"" , U=l , (A.IO) 

enjoys the following scaling symmetry: 

2 

(j) ^ cc/) , V ~^ V logc , (A. 11) 

a 

with trivial transformation rules for x, a, g, under which G — >■ c~^G, J — ^ J and the action remains 
invariant. The corresponding Noether charge is given by 

Q2 = a^e^/^GiT^W + - a^e-^gri . (A.12) 
a 

It is easy to check that eqs. (IA.7I) and (jA.12p are integrals of the equations of motion. Combining (jA.7p 
and (IA.12[) one finds that, in the model (jA.lOp . the following combination is also conserved 

\ a'^ ) a. 

By choosing the gauge a = r the equations (|A.3I) and (jA.4p become first order. Thus, in this gauge, one 
is left with a system of four coupled, independent, first order differential equations, (|A.3[) , (|A.4I) . (IA.7I) 
and (|XT3)) . 

One could try to use this first-order system to find, if not analytic, at least numerically simpler 
solutions for t] and 0. Unfortunately, this model does not seem to exhibit holographic superconductivity. 
The basic reason is that any solution asymptotic to AdS with t] ^ rjo (in particular, 77 — > 0) at infinity 
gives a mass to the photon; see eq. (|2.3|) . 



B AdS Reissner-Nordstrom solution via superpotential 

Here we provide a simple derivation of the AdS Reissner-Nordstrom solution using the superpotential 
method. We hope that this type of approach may also be useful for the search of analytic black hole 
solutions with non-trivial scalar hair r/. In general, one considers a Lagrangian of the type 

where Gy is a symmetric matrix, that may depend on the fields a' . The equations of motion are given 

by 



1/5 „ \da^ da'' dV 



d_ 

dr 

One then assumes that a superpotential W(a'') exists, which is such that it satisfies the following relation 



2\da'^^'") dr dr da' ' ^^''^^ 



1 dW dW 

If W is found, then one can show that the first order system, 

da'' dW 

automatically solves the equations of motion (the converse is in general not true; not all solutions of the 
second-order system solve the first-order system). 



13 



We now apply this procedure to our case. We consider the ansatz 

ds^ = -^dt^ + ^dr"^ + R^{dx^ + dy^) , A^dx" = 0(r) dt , (B.5) 
where S = S(r), R ~ R{r) and r/ = The Lagrangian then becomes 

V^C = 2R'S' + ^ G(77)0'' + U{ii) Usv'' + l-^J{v)c^' . (B.6) 



2 ^"^'12 '■"2 ' 2 S 
The equations of motion are given by 

R' S' ^ 1/2 3i?^^ 1 R? ^ ,o „N 

G0'2 + -r;' +_-[/ + -_ J</,2 ^ (B.7) 



i?5 4 5 ^ 4 ' 5- ' 4 52 

„ /2i?' d,-,Gi\ RJ 



^" + (f + f ) + 1 ^"^T^-s ^^^^ = « ■ (^-^^^ 

Combining eqs. p.l4p and p.lSp one can write the simpler relation 

_L i^R's)' = -1 - ^ G(,7)</''' • (B.ll) 

The relation (IB. 31) becomes 



9flTy5sVF + -^^^(9^14^)2 „ -^{d^Wf = l2R^U{i^) + J{i^)<f>^ . (B.12) 

A superpotential satisfying this relation can be found for 7/ = 0. We choose units such that L — 1. In 
this case (IB.12P reduces to (we assume U{0) — G(0) = 1) 

dnWdsW + -^{d^W)^ = 12R^ . (B.13) 
We have found two simple solutions to eq. (IB.13P 

W^=2{R^ + S)-^{<l>-f,)\ (B.14) 



W2 = - M^S- i - Ai)2 . (B.15) 

It is easy to see that the AdS Reissner-Nordstrom solution arises from these superpotentials. Consider, 
for example, the case of Wi. The first order system (jB.4[) then becomes 



R' = l, S' = 3R' - ^{^ - fif , 0' = -l(0-/,). (B.16) 



This is easily solved giving 

2 

R = r, S = r^-M + ^, (/> = ^ - E ^ (B.17) 

4r r 

where M and p are an integration constants. This reproduces the RNAdS solution (|2.26p . Clearly, it 
would be interesting to find a generalization of Wi or W2 to include 77 dependence. This could allow one 
to have more analytic control over the hairy black hole solution and thus over the properties of the phase 
transition. 
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